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We present a revised version of the so-called "yukawaon model", which was proposed 
for the purpose of a unified description of the lepton mixing matrix Upmns and the quark 
mixing matrix Vckm ■ It is assumed from a phenomenological point of view that the neutrino 
Dirac mass matrix Mo is given with a somewhat different structure from the charged lepton 
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_ mass matrix M e , although Mp> = M e was assumed in the previous model. As a result, the 

revised model predicts a reasonable value sin 2 2#i3 ~ 0.07 with keeping successful results for 
other parameters in Upmns as well as Vckm and quark and lepton mass ratios. 

> 
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CN " 1 Introduction 

In a series of papers [B El [3], the authors have investigated a unified description of the 
<3 lepton mixing matrix [1] Upmns and the quark mixing matrix [5] Vckm- The essential idea 

is as follows: (i) The Yukawa coupling constants Yt (/ = u, d, e, and so on) in the standard 
model are effectively given by vacuum expectation values (VEVs) of scalars ("yukawaon") Yy 
with 3x3 components, i.e. by (Yf)/A. Here A is an energy scale of the effective theory. 
(The yukawaon model is a kind of the "flavon" model [6j.) (ii) The model does not contain 
any coefficients which are dependent on the family numbers. The hierarchical structures of 
the effective Yukawa coupling constants originate only in a fundamental VEV matrix (^o)) 
whose hierarchical structure is ad hoc assumed and whose VEV values are fixed by the observed 
charged lepton masses, (hi) Relations among those VEV matrices are obtained from SUSY 
vacuum conditions for a given superpotential under family symmetries and R charges assumed. 
(Since we use the observed charged lepton mass values as the input values, it is a characteristic 
in the yukawaon model that adjustable parameters are quite few.) 

In the previous model[TJ[3], the quark and lepton mass matrices (charged lepton mass matrix 
M e , Dirac neutrino mass matrix Mo, down-quark mass matrix Mj, neutrino mass matrix M u , 
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and right-handed Majorana neutrino mass matrix Mr) are given as follows: 



M e = k e <S> (l + a e X 3 )$ , 
M D = M e , 

M d = k d [* (1 + a d X 3 )<Z> + m° d l] , 

M u = k' u M u M u , (1.1) 

M u = k u $ (l + a u X 3 )$ , 

My = M D M~ l Ml, 

M R = k R (M u M e + M e M u ) + ■■■ , 



where M e , $o, X 3 , ■ ■ ■ are 3x3 numerical matrices which results from VEV matrices of scalar 
fields. Here the VEV matrices <I>o> -^3> an d 1 have structures given by 
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The coefficients a/ (/ = e,u,d) which are important parameters in the model play an essential 
role in the mass ratios and mixings, while the family- number independent coefficients kf and 
k' u do not. The values of (x\,X2,x 3 ) with x\ + x\ + x\ = 1 are fixed by the observed charged 
lepton mass values under the given value of a e . (In an earlier model [7]j the charged lepton mass 
matrix M e was given by M e = k' e & e $ e and M d and M u are given by those with the replacement 
~~ ^e- The structures with (1 + afX 3 ) were suggested in a phenomenological model by 
Fusaoka and one of the authors [S].) 

The model has given almost successful unified description and predictions of Upmns and 
Vqkm- However, the model has failed to give the observed large mixing of 6*13 in Upmns'- the 
observed value is sin 2 2^13 ~ 0.09 [SJ, while the yukawaon model predicts sin 2 2^13 ~ 10~ 4 pQ. 
Since the model does not contain enough number of adjustable parameters as it is, it is hard to 
improve the prediction of sin 2 2#i3 ~ 10 -4 without the cost of other successful predictions. 

An introducing of the structure X2 into the yukawaon model has been done in Ref . [2] . In 
Ref.[2], the structure X2 is introduced in M e together with asummption Mp> = M e , and only the 
Upmns a r e discussed, but the predicted value of sin 2 26*13 is still small: sin 2 2^13 ~ 10~ 2 . The 
Vckm was not discussed in Ref. [2] . Even in a recent revised yukawaon model [3] , the predicted 
value was, at most, sin 2 26*13 ~ 0.03. 

In the present model, it is assumed from a phenomenological point of view that the neutrino 
Dirac mass matrix Mp> is given with a somewhat different structure from the charged lepton 
mass matrix M e : 

M D = k e $ (l + a D X 2 )^o, (1.3) 
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where 



Xo 




(1.4) 



Using this form we shall discuss Upmns as well as Vckm of the model. As to the structure X2, 
we will discuss in Sec. 2. When once we accept the form (1.3), we predict a reasonable value of 
sin 2 2^13 ~ 0.07 together with reasonable other parameters of Upmns, Vckm and quark and 
lepton mass ratios. 

2 Model 

We assume that a would-be Yukawa interaction is given as follows: 

Wy = ^e^HjH, + \v^t 3 R u + A^Y>| + ^d^qjH, + ^u^ qj H u , (2.1) 

where I = (vl-,&l) and q = (ul,cIl) are SU(2)i doublets. Under this definition of Yf, the 
CKM mixing matrix and the PMNS mixing matrix are given by Vckm = ulu d and Upmns = 
UeU u , respectively, where Uf are defined by ujM^MfUf = D 2 (Df are diagonal). In order to 
distinguish each yukawaon from others, we assume that Yf have different R charges from each 
other together with R charge conservation. (Of course, the R charge conservation is broken at 
the energy scale A'.) 

We assume the following superpotential for yukawaons by introducing fields 0* , P, E, E, 

E', E", E", E', cf> e , and d : 

W e = \ e + j(*o) ia (V'W + a e ±X ak E kl X^ 0^, (2.2) 

Wd = X {( E ')? Y D(E'f 3 + (*Zr (Eij + a D ±xT(E>y s X s ^ ($ )^} eg,, (2.3) 



w u = ^ {p^p^ + Y%E kl Y%} ei\ 



(2.4) 



W' = -± 
u A 



pkyupj + { Q o) ia (^ E n u + au J_x ak E kl X^ (<&£)«} Q% (2.5) 



i*oT ( (E"U + a d ^X ak E^X( p ) + m° d {E'UE"T P {E^ 



Q d . 



(2.6) 
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w R = | Mfl Yg + ^ [y^y^ + ^Y^yy + ^y^y^] } (2-7) 

Here, we have assumed family symmetries U(3)xU(3)'. The fundamental yukawaon $o is as- 
signed to (3, 3) of U(3)x U(3)', although quarks and leptons are still assigned to (3, 1) and 
yukawaons Yf are assigned to (6*, 1) of U(3)x U(3)'. In order to distinguish R charges between 
y e and Yd, we have introduced U(3)xU(3)' singlet scalar fields 4> e and 4> d - The present model is 
not anomaly free for U(3)xU(3)' if we consider only the fields which are given in Eqs.(2.1)-(2.7). 
In order to make the model anomaly free, we need further fields. However, since roles of such 
additional fields in the present model are, at present, not clear, we do not discuss such fields. 
From Eqs.(2.2) and (2.3) [and also (2.5) and (2.6)], we obtain 

R(E) + R(E) = R(E") + R(E"). (2.8) 

The VEVs of the introduced fields E, E, P, and P are described by the following superpotential 
by assuming R(EE) = R(PP) = 1: 

W E ,p = ^-Tt[EEPP] + ^-Tr[EE]Ti[PP], (2.9) 

which leads to 

(£}<£} ocl, (P)(P) ocl. (2.10) 
We assume specific solutions of Eq.(2.10): 

± {E ) = -(E) = 1, (2.11) 

VE VE 

— (P) = ^(P) j = diag(e-^,e-^, 1), (2.12) 
vp vp 

as the explicit forms of (E), (E), and (P). We assume similar superpotential forms for E" and 
E", and for E' and E'. 

From SUSY vacuum conditions dW/dQ = 0, we obtain the following relations: 

(Y e > = fc e (<D ) (1 + a e XX T ) (d#), (2.13) 

(Y D ) = k D (tf) (1 + a D X T X) ($ >, (2.14) 

(P)(Y u )(P) = k' u (Y u )(Y u ), (2.15) 

(Y u ) = fc u (cD ) (1 + a u XX T ) ($1), (2.16) 

(Y d ) = k d [($ > (1 + a d XX T ) (tf) + m° d l] , (2.17) 
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(Y R ) = k R ((Y e )(Y u ) + (Y u )(Y e )+£(Y D )(Y D )) , (2.18) 

where, for convenience, we have already put (E) as 1, and so on. Here, since we have assumed 
that all fields take (0) = 0, we do not need to consider vacuum conditions for other fields 
dW/dY e = 0, because those alway contain (9). Thus, mass matrices are given by M e = (Y e ), 
M D = {Y D }, M u = k' u M u M u , M u = (Y u ), M d = (Y d ), M v = M D M R X M T D , and Mr = (Y R ). 
The most curious assumption is to assume the VEV matrix form of the scalar X as 



1 , , 1 
— X = - 



/ 1 1 \ 

1 1 
V 1 1 o J 



(2.19) 



At present, there is no idea for this form. This is only ad hoc assumption. However, as seen 
later, we can obtain good fitting for the neutrino mixing angle sin 2 29\ 3 due to this assumption. 
The form (2.19) leads to 

(X){X T ) = \x 3 , (X T )(X) = ^X 2 , (2.20) 

together with (X)(X) = (X), where X% and X 2 is defined by Eqs. (1.2) and (1.4), respectively, 
and, for simplicity, we have put vx = 1 because we are interested only in the relative ratios 
among the family components. 

3 Parameter fitting 

We summarize our mass matrices in the present model as follows: 



M e = k e $ (l + a e X 3 )$T , (3.1) 

M D = k D ^(l + a D X 2 )^ , (3.2) 

PM U P = k' u M u M u , (3.3) 

M u = k u $> (1 + a u e ia "X 3 ) (3.4) 

M d = k d [$ (1 + a d X 3 )^ + m°l] , (3.5) 

M u = M D M R l Ml, (3.6) 

M R = k R [M e M u + M u M e + £m d Md) , (3.7) 



where, for convenience, we have dropped the notations "(" and "}". Since we are interested 
only in the mass ratios and mixings, hereafter, we will use dimensionless expressions $o = 
diag(xi, x 2 , x 3 ), P = diag(e l( ^ 1 , e 1 ^ 2 , 1), and E = diag(l, 1, 1). For simplicity, we have regarded 
the parameter a d as real correspondingly to the parameter a e . The parameters are re-refined by 
Eqs.(3.1)-(3.5). 
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In the present model, we have 9 adjustable parameters [a e , cld, (a 
and for the 16 observable quantities (6 mass ratios in the up-quark-, down-quark-, and 
neutrino-sectors, 4 CKM mixing parameters, and 4+2 PMNS mixing parameters). 

Since we do not change the mass matrix structures for M e , M u , and from the previous 
paper [3], we use the following parameter values of a e and (a u ,a u ) 

(a e ,a u ,a u ) ~ (8, -1.35, ±8°), (3.8) 

which are fixed from the observed values of m c /mt, m u /m c , and sin 2 26*23: 

rh = ^ = 0.0451°;™, r& = ^ = 0.06 ± 0.005, (3.9) 

at {i = mz |10| . and sin 2 26*23 > 0.95 |11| . (These values will be fine-tuned in whole parameter 
fitting of Upmns and Vckm later.) Note that the neutrino sector of the model is different from 
the previous model, however the predictrd values of sin 2 26*23 are almost the same before. 
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Figure 1: Lepton mixing parameters sin 2 26*i2, sin 2 26*23, sin 2 26*i3, and the neutrino mass squared 
ratio R v versus the parameter cld- ("solar", "atm", "tl3", and "10 R" denote curves of sin 2 29x2, 
sin 2 2023, sin 2 26*13, and R u x 10, respectively. Other parameter values are taken as a e = 7.5, 
a u = —1.35, and a u = 7.6°. 

First, let us investigate lepton sector. In the revised model, a new parameter cld is added. 
We illustrate the behaviors of Lepton mixing parameters sin 2 26*12, sin 2 26*23, sin 2 26*i3, and the 
neutrino mass squared ratio R u versus the parameter for a case of £° = 0. As seen in Fig.l, 
the parameter cld does not change the prediction sin 2 26*23 ~ 1 in the previous model. Also, note 
that the prediction of sin 2 26*i3 is insensitive to the parameter cld, i-e. sin 2 26*13 ~ 0.08. Only 
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the predictions of sin 2 26*i2 and R v = {m\ ~ m i)/( m ! ~~ m |) are sensitive to the parameter op. 
In order to fit the observed value p] sin 2 26 12 = 0.857 ± 0.024, we take ao = 9.01. However, in 
the model with £° = 0, the value ao = 9.01 cannot fit the observed value [11] of R v , 

K= H _ (7-50 ± 0.20) x 10-5 eV 2 _ . 14 

^-Aml 2 " (2.32l^ 2 ) x 10-3 eV 2 " ^ M -aw) x 10 ■ ^ 10 J 

The non-zero parameter has phenomenologically been brought in order to adjust the predicted 
value of R u . The predicted values of sin 2 26*23, sin 2012, and sin 2 26\s are almost unchanged 
against the parameter In order to fit the neutrino mass ratio R u , we take £°, = —0.78. 

Next, we discuss quark sector. Since we have fixed the five parameters a e , a u , a u , op, 
and t®, we have remaining four parameters for six observables (2 down-quark mass ratios and 4 
CKM mixing parameters). The parameters and m d are used to fit the observed down-quark 
mass ratios [10] 



d ms _ +0.006 d = _ n nC) „ + o.005 /o i i \ 

23 _ ~ -0-006' r 12 - — - u - UD,; >-0.003' V^- 11 ; 

respectively. Therefore, the four CKM ming parameters are described only by two parameters 
(01 5 02)- As shown in Fig. 2, all the experimental constraints on CKM mixing matrix elements 
are satisfied by fine tuning with use of only two parameters (0i, 4*2) ■ 

Finally, we do fine-tuning of whole parameter values in order to give more improved fitting 
with the whole data. Our final result is as follows: under the parameter values 

a e = 7.5, ao = 9.01, (a u , a u ) = (-1.35, -7.6°), a d = 25, 

m° d = 0.0115, C° = -0.78, (</)i,0 2 ) = (177.0°, 197.4°), (3.12) 



we obtain 



= 0.0465, r& = 0.0614, rf 2 = 0.0569, = 0.0240, (3.13) 



sin 2 26*23 = 0.969, sin 2 26 12 = 0.860, sin 2 20 13 = 0.0711, R v = 0.0324, (3.14) 
b l cp = -131° {f = -2.3 x 10~ 2 ), (3.15) 

\(V C km)us\ = 0.2271, \{V C KM)cb\ = 0.0394, \{V CK M)ub\ = 0.00347, |(Vcxm)wI = 0.00780, 

(3.16) 

S q cp = 59.6° (J q = 2.6 x 10~ 5 ). (3.17) 

Here, 5^ CP and 5^,p are Dirac CP violating phases in the standard conventions of Upmns and 
Vckm, respectively. 

The predicted value of sin 2 26*13 — 0.0711 is still somewhat small compared with the observed 
value sin 2 26*13 = 0.098 ± 0.013 [IT]. So far, we have assumed that the parameter £j) is real. If 
we consider that the parameter £° is complex, £® — > £,®e ta ", we can adjust the value of sin 2 26*13 
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Figure 2: Contour plots in the (<j)i, 02 ) parameter plane, which are shown by using experimental 
constraints on CKM mixing matrix elements: \{Vckm)us\ = 0.2252 ± 0.0009, |(Vc , a:m)c6| = 
0.0406 ±0.0013, \{V CK M)ub\ = 0.00389 ± 0.00044, and \{V CK M)td\ = 0.0084 ± 0.0006. The CKM 
elements depends on only the parameter set of [a e , (a u , a u ), ad, m d , 0i, and 02]. Here we 
present contour plots of the CKM elements in the (0i, 02 ) parameter plane by taking the values 
of other parameters as a e = 7.5, a u = —1.35, a u = —7.6°, = 25, and m d = 0.0115. We find 
that (0i, 2 )=(177.O°, 197.4°) is consistent with all the CKM constraints. 

without changing other predicted values as seen in Fig. 3. However, such a modification by the 
parameter a v is not essential in the present model, so that we will regard that the parameter £® 
is real. 

We can also predict neutrino masses, for the parameters given (3.12) with real 

m ul ~ 0.0048 eV, m u2 ^ 0.0101 eV, m u3 ~ 0.0503 eV, (3.18) 

by using the input value [T2] Am| 2 ~ 0.00243 eV 2 . We also predict the effective Majorana 
neutrino mass [13] (m) in the neutrinoless double beta decay as 

(m) = \mx(UpMNs)li + m 2 {UpMNs)l2 + m 3 (UpMNs) 2 e 3 \ ~ 7-3 x 10~ 4 eV. (3.19) 

6 Concluding remarks 

In conclusion, by assuming VEV matrix forms of the yukawaons Eqs.(3.1) -(3.7), we have 
obtained reasonable Vckm and Upmns mixing parameters together with quark and neutrino 
mass ratios. The major change from the previous yukawaon models is in the form of Mp>. 
Although we give the form by assuming the VEV matrix X which is geven by Eq.(2.19), and 
by considering the mechanism X T X versus XX T , it is still phenomenological and somewhat 
factitious. However, when once we accept the form of Md, we can obtain sin 2 2#i3 ~ 0.07 whose 
value is not sensitive to the other parameters. 
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Figure 3: Lepton mixing parameters sin 2 26*i2, sin 2 26*23, sin 2 26*i3, and the neutrino mass squared 
ratio R v versus the phase parameter a v . ("solar", "atm", "10 tl3", and "10 R" denote curves of 
sin 2 20i2, sin 2 2023, sin 2 26*i3 x 10, and R v x 10, respectively. Here the a u dependence is presented 
under the parameter values given by (3.12). 

Note that the present model does not have any family-dependent parameters except for 
(xi, X2, X3) in (<J>o) and {<j>\, 02 ) in (P)- The parameter values (xi, X2, £3) have been fixed by the 
observed charged lepton masses. Therefore, the model have only 9 adjustable parameters for 16 
observables. The 5 parameter values of 9 parameters, (a e , an, (a u , a u ), and have been fixed 
by the observed values m u /m c , m c /m t , sin 2 26> 12 , sin 2 26*23, an d R v . Especially, the parameter 
an has been fixed the observed value sin 2 20 12. The parameter £® has been introduced only in 
order to adjust the ratio R v . (In other words, the value of £® has been fixed by R° bs , the value 
(3.10).) Logically speaking, we need four observed values in order to fix the four parameters 
a e , od, and (a u ,a u ). However, as seen in Fig.l, the values sin 2 2023 ~ 0-9 and sin 2 26*13 ~ 0.08 
are almost determined independently of the parameter ao when we fix a e and (a u , a u ) from the 
observed up-quark mass ratios. Therefore, sin 2 26*23 ~ 0-9 and sin 2 26*13 ~ 0.08 can substantially 
be regarded as predictions in the present model. Of course, after we fix the 5 parameters, 
predictions are 6 quantities: sin 2 26*13, 2 neutrino mass ratios, CP-violating phase parameter, 
and 2 Majorana phases. 

Of the remaining 4 parameters ad, rn d , and (01,02), the parameters ad and m d are deter- 
mined by the down-quark mass ratios m s /mb and m,d/m s , respectively. Therefore, the 4 CKM 
mixing parameters are predicted only by adjusting the two parameters ((pi, <p2)- We can obtain 
reasonable predictions of the CKM mixing parameters. 

The present model is still in a level of a phenomenological model. Nevertheless, it seems 
that the yukawaon model offers us a promising hint for a unified mass matrix model for quarks 
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and leptons, because quark and lepton mass matrices are given each other by their mixed terms. 
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